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ABSTRACT 


A system  of  linear  electroelastic  equations  for  small  fields  superposed 
on  a bias  is  applied  in  the  determination  of  the  velocity  of  acoustic  surface 
waves  in  piezoelectric  substrates  subject  to  flexural  biasing  stresses.  The 
influence  of  the  biasing  stresses  appears  in  the  boundary  conditions  as  well 
as  the  differential  equations.  Direct  calculations  performed  for  both  quartz 
and  lithium  niobate  when  the  spatial  variation  of  the  flexural  biasing  state 
is  omitted  indicate  that  the  biasing  stresses  in  the  boundary  conditions  have 
an  important  influence  on  the  surface  wave  velocity.  In  addition,  perturba- 
tion calculations  are  performed  which  include  the  influence  of  the  spatial 
variation  of  all  flexural  biasing  terms  and  it  is  shown  that,  for  substrate 
thickness-to-wavelength  ratios  well  within  the  practical  range,  the  spatial 
variation  in  the  biasing  state  has  an  appreciable  effect  on  the  velocity  of 


acoustic  surface  waves. 


1.  Introduction 


In  recent  years,  a number  of  workers  have  measured  the  change  in 

acoustic . surface  wave  velocity  due  to  applied  static  biasing  stresses, 

2 

primarily  for  application  as  pressure  sensors.  In  one  of  the  aforementioned 

works  an  analytical  treatment  was  presented,  whose  validity  seems  questionable 

to  us  because  the  treatment  is  not  related  to  the  fundamental  equations. 

4 

More  recently,  another  independent  analytical  treatment  was  presented  , which 

properly  considers  the  elastic  nonlinearities  in  the  differential  equations 

but  not  in  the  boundary  conditions. 

In  this  paper  a system  of  linear  electroelastic  equations  for  small 

5 

fields  superposed  on  a bias  , which  were  obtained  from  general  nonlinear 
rotationally  invariant  electroelastic  equations^,  is  applied  in  the  determina- 
tion of  the  velocity  of  surface  waves  in  piezoelectric  substrates  subject  to 
flexural  biasing  stresses.  The  influence  of  the  biasing  stresses  appears  in 
the  boundary  conditions  as  well  as  the  differential  equations.  Although  the 
electric,  electroelastic  and  elastic  nonlinearities  are  included  in  the  basic 
equations,  only  the  elastic  nonlinearities  are  included  in  the  calculations 
because  the  nonlinear  electroelastic  coefficients  are  not  known  for  either 
lithium  niobate  or  quartz,  the  two  materials  considered.  The  nonlinear 
electroelastic  coefficients  may  very  well  be  significant  in  lithium  niobate, 
even  though  they  are  probably  negligible  in  quartz.  Nevertheless,  numerical 
results  are  obtained  for  Y-Z  lithium  niobate  and  y-cut  and  ST-cut  quartz  using 

the  published  values  of  the  second  and  third-order  elastic,  piezoelectric  and 

7 8 

dielectric  constants  for  these  materials  ’ . 

The  special  case  of  biasing  stresses  due  to  flexural  loading  of  a simply 
supported  plate  is  considered  in  detail  and  all  components  of  the  static  deforma- 
tion that  occur  in  the  linear  equations  for  small  fields  superposed  on  a 


2. 


bias  are  analytically  related  to  the  applied  loading.  Analyses  are  performed 
for  both  rotated  Y-cut  quartz  and  lithium  niobate  subject  to  biasing  stresses 
both  in  and  normal  to  the  propagation  direction.  Numerical  results  are  pre- 
sented for  both  Y-cut  and  ST-cut  quartz  for  propagation  in  the  digonal  direc- 
tion and  for  Y-cut  lithium  niobate  for  propagation  in  the  Z-direction. 

9 

Finally,  a perturbation  equation  for  the  change  in  surface  wave  velocity 
under  any  biasing  stress  is  employed  and  shown  to  yield  results  indistinguishable 
from  those  obtained  from  the  direct  calculation  for  the  particular  case  of 
flexure  considered  here  for  wavelengths  very  small  compared  to  the  thickness 
of  the  plate.  In  fact,  since  the  perturbation  procedure  can  readily  treat 
spatially  varying  biasing  states,  for  which  a direct  calculation  cannot  be 
performed,  it  has  a significant  advantage  over  a direct  calculation  in  the 
determination  of  the  small  changes  in  surface  wave  velocities  due  to  biasing 
states.  Indeed,  in  the  case  of  flexural  loading  considered  here  the  biasing 
terms  actually  are  functions  of  position  and  this  fact  is  ignored  in  the  direct 
calculation.  When  spatial  variation  of  the  biasing  terms  is  omitted  in  the 
perturbation  calculation,  the  aforementioned  agreement  with  the  direct  calcula- 
tion for  the  thick  substrate  results.  However,  when  the  spatial  variation  in 
the  biasing  terms  is  included  in  the  perturbation  calculation,  a readily  dis- 
tinguishable change  from  the  results  of  the  direct  calculation  occurs  for  substrate 
thicknesses  less  than  about  100  wavelengths.  The  deviation  increases  signifi- 
cantly with  decreasing  thickness  and  becomes  appreciable  for  thicknesses  of  the 
order  of,  say,  50  wavelengths,  which  is  well  within  the  practical  range.  The 
results  of  the  calculation  for  ST-cut  quartz  subject  to  a biasing  stress  in  the 
propagation  direction  have  been  compared  with  experiment10  and  the  agreement  is 
quite  good. 
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3. 


2.  Linear  Electroelastic  Equations  for Small  Fields  Superposed  on  a Bias 


The  linear  electroelastic  equations  for  small  fields  superposed  on  a bias 


may  be  written  in  the  form 


where 


\.Y,L  = P°V  \l=°’ 

\y  = + 2MLY^,M’ 

\ ~ RLMQ'U0,  m + rlm^,  m ’ 


G = C + C , R — — G — 6 . 

ILNMV  2L"NMv  2LYMv’  2LM  LM  LM  ’ 

2MLY  = iMLY  = eMLY  + SMLY  * 

2LYMv “ TLM6YV  + 3LYMvABEAB  + 2LYKMWv,  K 

-1  1 

+ c w — ic  CP+a 

2LKMV  Y,  K iALYMv  ,A  yLYM\)  ’ 


1 ~1  2 
E -e  w -b  cp+q 
IC  BC  MLK  Y,  K AMLY  ,A  yMLY  * 


*•  1 -1  11 

e = b E - X_  cp  - 2e  J E , (2.3) 

LM  LMCD  CD  £LMCy,C  O ML  ’ ' 
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and  g_  and  g are  lengthy  expressions  ’ involving  the  biasing  electric 
L tMv  ML  Y 

field  and  static  deformation,  which  vanish  when  the  biasing  electric  field 

vanishes  and  are  not  of  interest  here.  Consequently,  we  do  not  present  them 

here.  Equations  (2.1)  constitute  the  stress  equations  of  motion  and  charge 

equation  of  electrostatics  referred  to  the  position  coordinates  of  material 

points  before  the  static  deformation  is  applied,  which  are  called  reference 

coordinates  and  are  denoted  X . Equations  (2.2)  are  the  linear  electroelastic 

M 

constitutive  equations  and  Eqs.  (2.3)  contain  the  definitions  of  the  effective 
coefficients  defined  therein.  In  (2.1)  - (2.3)  K^,  u^  and  2$^  denote  the  com- 
ponents of  the  small  field  Piola-Kirchhoff  stress  tensor  which  is  asymmetric, 

5 

mechanical  displacement  vector  , and  reference  electric  displacement  vector  , 
respectively;  p°  and  9 denote  the  reference  mass  density  and  small  field 


4. 


electric  potential,  respectively;  c , e and  e denote  the  second  order 

2 L yMV  ML  t LM 


elastic,  piezoelectric  and  dielectric  constants,  respectively,  which  are  the 
constants  that  occur  in  the  ordinary  linear  theory  of  piezoelectricity. 

Before  discussing  the  remaining  quantities  in  (2.3)  we  note  that  in  a 
motion  the  present  position  of  material  points  y may  be  written 


"2+X<V  +S(XL>t)> 


(2.4) 


where  w denotes  the  displacement  due  to  the  applied  static  loading.  In  addi- 
tion, we  note  that  the  electric  potentials  $ in  the  body  and  in  free  space 
may  be  written  in  the  respective  forms 


$ (XL,  t)  = 91  (XL)  + (XL,  t)  , ♦ (XL,  t)  = t1  (XL)  + * (XL,  t)  , 


(2.5) 


~1  A1 

where  cp  and  f denote  the  static  biasing  electric  potentials  in  the  body  and 


in  free  space.  In  (2.3)  and  E?_  denote  the  components  of  the  static 

LM  AB 


biasing  stress  and  strain,  respectively,  and  J is  the  Jacobian  of  the  static  de- 
formation. The  biasing  variables  satisfy  the  appropriate  static  equations 
given  in  (66)  - (72)  of  Ref. 5,  or  the  equivalent  equations  using  reference 
coordinates  as  independent  variables.  Since  we  are  interested  in  small  biasing 
strains  only,  we  have 

_1  1 . .1 

(2.6) 


EAB  (WA,B  +WB,A}>  jl~1> 


TLM  ?LMRSWR,  S + eRLMP,  M * 


(2.7) 


The  quantities  b^,  and  denote  the  third-order  elastic, 

electrostrictive,  third-order  electric  permeability  and  first-order  electro- 


elastic constants,  respectively,  and  Gq  denotes  the  electric  permittivity 


of  free  space.  For  obvious  reasons  the  notation  employed  here  is  designed  to 
be  consistent  with  the  notation  of  Refs. 5 and  9.  The  carets  over  many  variables 
have  been  employed  here  because  we  consistently  use  the  reference  coordinates 


5. 


as  independent  variables.  In  particular  the  Greek  subscripts  are  employed  in 
order  to  be  consistent  with  the  notation  of  Refs. 5 and  9,  in  which  Greek  sub- 
scripts were  consistently  employed  to  denote  the  intermediate  coordinates  £ . 
However,  the  Greek  subscripts  are  not  really  important  for  any  of  the  opera- 
tions performed  here  and  may  be  used  interchangeably  with  the  capital  Latin 
subscripts  which  denote  reference  coordinates  X . The  cycles  above  variables 
also  have  been  introduced  for  consistency  with  Refs. 5 and  9.  We  have  employed 
Cartesian  tensor  notation  and  introduced  vector  notation  in  Eq.  (2.4)  only, 
and  the  convention  that  a comma  followed  by  an  index  denotes  partial  differ- 
entiation with  respect  to  a reference  coordinate,  the  dot  notation  for  partial 
differentiation  with  respect  to  time  and  the  summation  convention  for  repeated 
tensor  indices. 

. 5 g 

To  these  equations  we  must  adjoin  the  boundary  conditions  ’ 

= 0> 

<P  = VLUL  + I LWL  + t (2.8) 

where  w and  u in  (2.8),,  respectively,  denote  the  static  biasing  and  small 
field  dynamic  components  of  the  mechanical  displacement  at  the  surface  of  the 
solid,  N denotes  the  unit  normal  to  the  reference  position  of  the  surface  and 

~f  f f f ? 

L,y  ~ ®LY3Up  + fL'VM0U0,M  + gLYM  , M ’ 

^•^US*lLv"v,M-eo*,L  (2'9> 

The  coefficients  G^,  G*w  G^,  R^  and  R^  are  lengthy  expressions5’ 9 
involving  the  biasing  electric  field  and  static  deformation,  which  vanish  when 

the  biasing  electric  field  vanishes  and  are  not  of  interest  here.  Consequently, 

~ 59 

we  do  not  present  them  here.  In  free  space  i|i  satisfies  Laplaces  equation  ’ , 

i.e. , 


6. 


it  = 0 . 
\ KK 


(2.10) 


Since  the  nonlinear  electroelastic  coefficients  are  not  known  for  either 
lithium  niobate  or  quartz,  or  to  our  knowledge  for  any  material,  only  the 
elastic  nonlinearities  will  be  included  in  the  remainder  of  the  analysis  and 
in  the  calculations.  However,  we  wish  to  emphasize  the  fact  that  we  could 
equally  readily  have  included  the  electroelastic  nonlinearities  if  the  coef- 
ficients were  known.  In  the  case  of  a purely  elastic  bias  Eqs.  (2.1)  and  (2.2) 


remain  unchanged,  Eq.  (2.3)^  reduces  to  the  form 

2LYM\>  ~ TLM6YV  + 3LVMVABEAB  + 2h\KMWy,  K + 2LKMVWY,  K 


and  Eqs.  (2.3)  and  (2.3).  are  replaced  by 
5 o 


eMLY  ° ’ ®LM  ° * 


(2.11) 


(2.12) 


In  the  case  of  the  purely  elastic  bias  the  boundary  conditions  in  (2.8)  reduce  to 

VlV  - 0 . SlA  ■ - Vo*,  !>*-♦>  (2- 131 


and  the  free  space  equation  (2.10)  remains  unchanged. 


3.  Static  Flexural  Loading 

Since  an  acoustic  surface  wave  is  sharply  confined  to  the  vicinity  of  the 

. 

surface  of  the  substrate  in  which  it  propagates  and  in  static  flexure  of  a thin 
plate  the  magnitude  of  the  stress  is  maximum  at  the  outer  fibers,  it  is  ad- 
vantageous to  employ  flexural  loading  structures  in  order  to  maximize  the  in- 
fluence on  the  surface  wave  velocity  for  a given  loading.  Consequently,  in 
this  section  we  briefly  discuss  cylindrical  flexure  of  thin  plates  for  the  type 
of  loading  shown  in  Fig.l.  The  elementary  theory  of  cylindrical  flexure  of 


7 


thin  plates,  which  is  known  to  be  extremely  accurate  when  applicable,  rests 
on  a few  critical  simplifying  assumptions.  In  most  treatments  of  elementary 
flexure,  the  theory  is  not  presented  in  a way  that  readily  indicates  the 
relation  between  the  flexural  variables  employed  and  the  mechanical  displace- 
ment vector  of  the  three-dimensional  t.ieory,  which  is  needed  here  for  the 
static  biasing  terms  in  Eq.  (2.11),  nor  is  the  theory  discussed  in  the  case 
of  anisotropy.  However,  Mindlin^1  has  presented  a derivation  of  the  equations 
of  the  theory  of  flexure  of  thin  anisotropic  plates  in  a manner  that  clearly 


indicates  the  relation  between  the  plate  variables  and  the  mechanical  dis- 
placement vector  of  the  three-dimensional  equations.  The  brief  discussion 
presented  for  the  special  case  of  interest  here  is  based  on  material  in 
Ref. 11. 


In  the  case  of  cylindrical  flexure  in  the  x^-direction  with  restrained 
motion  in  the  x^-direction  the  aforementioned  critical  assumptions  are  that 
the  static  displacement  field  w^  may  be  written  in  the  form 

= z2wi1}  (zi}  , w2  = w20)  (zi)  + Z2W2(2)  (Z1}  » w3  = 0>  (3.1) 

11  (o) 

the  zero  order  shearing  strain  vanishes  while  the  zero  order  shearing 

stress  T^°^  is  related  to  the  bending  moment  T^^  = M by  the  usual  static 


equation 


T(0)  = M 
12  , 1 » 


and  the  first  order  thickness  stress  vanishes.  The  conditions  = 0 and 

T22^  = respectively,  yield  the  important  relations 

“a2’-  <cl/2c2X"  ' <3'3’ 


Substituting  from  (3.3)  into  the  constitutive  relation  for  M,  we  obtain 


'2h3/3>cn“2°u 


8 


* 2 . 
cn  cn  ” ci  > 


(3.5) 


where 

''ll  - "ll  ~ v'12/v'22 
and  it  is  to  be  emphasized  that  these  equations  hold  for  arbitrary  orienta- 
tions of  in  the  anisotropic  case.  From  the  loading  shown  in  Fig.l,  in  the 
region  in  which  the  surface  wave  velocity  is  measured,  i.e.,  \z^\  < b,  we  have 


M=  - Pa 


(3.6) 


where  P is  the  load  per  linear  dimension,  which  with  (3.2)  indicates  that 
= 0.  Since  the  slope  of  the  deflection  curve  vanishes  at  = 0,  from 
(3.4)  we  obtain 

(o)  3 

(3.7) 


w2(°|  = - zx  (3/2h3)M/c*1  , 


which  relates  the  loading  to  the  slope  of  the  deflection  curve.  From  (3.1) 
and  (3.3)  we  obtain  the  important  relation 


W2,2  = “ (C12/C22)W1,1' 


Now,  from  (3.1),  (3.3)  and  (3.7),  we  obtain 

W1  1 = z2(3/2h3)M/ci!  • 


Moreover,  from  (3.1),  (3.3)  and  (3.7),  we  find 

3M  3M 


w,  „ = z 


W„  = - z. 


1.2  1 „ 3 * > 2.1  1 3 


(3.8) 


(3.9) 


(3.10) 


'11  ““  ~11 
and  all  other  three-dimensional  displacement  gradients  vanish.  From  (2.6), 
(3. 1)  ^ and  (3.7)-  (3.10),  we  obtain 

E”  = z~  ~ > e22  = _ z2  — ^rrir1  > E^  = 0>  E^v=°>  f3-11* 


11  2 ov, 

2h  cu 


2 2h3c^1  C22  ' “12 


3K 


for  the  three-dimensional  strain  field.  In  addition,  from  the  usual  three- 
dimensional  linear  anisotropic  stress-strain  relations  and  (3.11),  we  have 


1 3M  1 

T 3 z — ■ T 

2 3 » A33 

2n 


_1  - Tl1  f 

» T33  * V 


C13  " ' 


°23C12 


N 1 T11  ( 

/ » T31  ” _*  '15  ' 


C25C12> 


(3.12). 


which  are  the  nonvanishing  components  of  stress  resulting  from  (3.11). 

Note  that  at  the  top  fiber  = - h,  is  a maximum  and  has  the  value 

T11  = 3M/2h2  . (3.13) 

In  the  case  of  cylindrical  flexure  in  the  x^-direction,  the  analytical 
treatment  is  the  same  except  for  the  fact  that  the  index  3 replaces  1 through- 
out. Under  these  circumstances  the  pertinent  displacement  gradients  and 
strain  components  take  the  respective  forms 


B,,  = z- 


W3,  2 

Z3  -.3  * * 

2h  C33 

W2,3 

3M 

E1  ' = - z 

3M 

3 * 

’ E22  2 

„ 3 : 

2h  c 

33 

2h  c. 

E„  _ = 0 . E._  = 0 


(3.14) 


(3.15) 


in  place  of  the  forms  in  (3.10)  and  (3.11)  and  where 


C33  C33  “ C32/C22  * 


In  addition,  in  place  of  (3.12)  and  (3.13)  we  have,  respectively, 


T1  M 33  | 

< _ C12C32^ 

T1 

1 T1  = -H  1 

f _ C25C32\ 

„ 3 > ail  * ' 

2h  C33 

^ c22  J 

1 ’ ri3  ~ 1 
°33 

^35  c22  ) 

T*™  = 3M/2h2 


(3.16) 


(3.17) 


(3.18) 


Thus,  depending  on  the  loading  situation,  we  use  either  Eqs.  (3.10)  - (3.12)  or 
Eqs.  (3.14),  (3.15)  and  (3.17)  in  the  static  bias  equation  (2.11). 

4.  Application  to  Rotated  Y-Cut  Quartz  and  Y-Cut  Lithium  Niobate 


For  the  case  of  static  flexural  biasing  stresses  both  in  and  normal  to  the 
propagation  (X^) -direction,  for  rotated  Y-cut  quartz,  Eqs.  (2.1)  and  (2.2),  with 


(2.11),  the  appropriate  relations  in  Sec. 3 and 


W2,2a’“  (C32/C22)W3,3  ’ 


take  the  form12’ 13 


C11U1, 11 + {C12  + C66 )U2, 12  + (C14  + C56 )U3, 12  + C66U1, 22  + ®11^, 11  + ®26^, 22 
+ EU[(CU  + 2C11+C111  + QC112)U1,11+  (2c66+C166+aC266)ul,22 

+ (C12  +c66  +C112  +C166  +a(cl2  +c66  + C122  + C266}  )u2, 12  + (c14  + C56 
+ c114  + C156  + a (C124  + C256 ) >U3, 123  + Wl,  2 [c66U2, 11 + °22U2, 22  + 2 (C12 
+ C66)U1,12 +C56U3,ll'rC24U3,  223  + W2, 1 fCllU2, 11  + °66U2,  22J 
+ E33 [ (C13  + C113  + * (C12  + C112) )U1, 11 + (C123  + C366  + 3 (°12  + C122  + °66 

+ C266 3 )U2, 12  + (C14  + C56  + C134  + C356  + P (C124  + )U3, 12  + (c366+  . 

+ Pc266)Ul,22]+W3,2(C12  + C66)U3,12+W2,3(C14  + C56)u2,12  = p{il  ’ 

C56U3,11  + (C12  +C66)U1,12  +C66U2,11  + C22U2,22  + C24U3,  22  + (®26  + e 12^,12 
+ E11 C (c12  + c66  + C112  + C166  + a (C12  + c66  + C122  + °266 5 )U1, 12  + (cll  + °166 
+Q'(2C66+C266))U2,11+  (C122+Q,(2C22+C222))U2,22+  (c156+“(c56 
+ C256))U3,11+  (°124 +Q,(C24 +C224))U3,22]  + "l,  2 [°66U1, 11  + C22U1,  221 

+ W2, 1 [C11U1, 11 + °66U1, 22  + 2 (C12  + C66 3 U2, 12  + (C14  + C56 )u3, 121 
+ £33  [ (c13+c366  + 3 (c12  + 2c66  + c266)  )u2>  L1  + (c123  + c366  + 3 (c12  + c66  + c122 

+C266))U1,12  + (C56  +C356  + P (°56  + C256 ) )u3, 11  + (°223  + P (2c22  + C222) )U2, 22 

+ (C24  + C234  + P (C24  + C224) )u3, 221  + *3, 2 [C22U3, 22  + C66U3, ll1  + W2, 3 1 (C14 

+ C56)U1,12  + 2C56U2,11  + 2C24U2,22+C55U3,11  + C44U3,221  = P*2  ’ 


11 


C55U3, 11  + (C56  + C14)U1, 12  + C56U2, 11  + C24U2,  22  + °44U3,  22  + (625  + C14)CP,  12 


+ E11 [ (C14  + C56  + C114  + C156  + “ (C124  + C256  > > Ul, 12  + (C156  + “ (C56  + C256  > > U2, 11 


+ (C124  + “ (C24  + C224) )U2, 22  + (C11 + C155  + “C255)U3, 11 + (cl44  + ac244)u3, 221 


+ Wl,2tC56Ul,ll+C24Ul,22]  + W2, 1 (°14  + C56  >U2, 12  + E33  [ (C14  + C56  + °134 


+ C356  + P (C124  + °256 ) )ul, 12  + (°56  + C356  + P (°56  + C256 } )u2, 11 + (C24  + C234 


+ 3 (c24  + c224>  )u2> 22  + (C13  + Pc12  + 2°55  + C355  + ^C255)U3, 11 + (2c44  + C344 


+ PC244)U3,22]  +W3,2[(C12  + C66)U1,12+C66U2,11  + C22U2,22  + 2C56U3, 
+ 2C24U3,22]  +W2,3[C44U2,22  + C55U2,111  = 0*3  > 


11 


ellUl, 11 + (el2  + ® 26 )U2, 12  + (®14  + ®25)U3, 12  + ®26U1, 22 


where 


- 6llVl  - e22^,22  - ° > 

a = ~ C12/C22>  ^ = _C32/C22* 


(4.2) 


(4.3) 


From  (2.13),  with  (2.10),  for  straight-crested  waves  with  propagation  wave- 
number  |,  we  have  the  boundary  conditions 


K2y  = °,  ^2  = " eo^>  at  X,  = 0, 


(4.4) 


2 o T ’ — 2 

12  13 

and  from  (2.2)  and  (2.11)  we  obtain  ’ the  pertinent  constitutive  relations 


K21 = C56U3, 1 + C66  (U1,  2 + U2, 1>  + 6 26^,  2 + E11 [ (2c66  + C166  + ^266  )ul,  2 


+ (c66  +C166  +<y(c66  +C266))U2,1+  (c56  + c156  +ac256)u3,l]  +wl,2[(°12 


+ C66 )U1, 1 + C22U2, 2 + C24U3, 21  + *2, 1C66U2, 2 + E33 [ (C366  + Pc266 5 ul, 2 


+ (C366  +P(C66+C266))U2,1+  (c56  + C356  + Pc256 )u3, 1]  + W3, 2C66U3, 1 + W2,  3C56U2, 1 ’ 


r 


11. 


C55U3, 11  + (C56  + C14)U1, 12  + C56U2, 11  + °24U2,  22  + °44U3,  22  + (®25  + ei4)Cf>,  12 

+ E ^ [ <c14  + c56  + c1M  + c156  + a (c124  + c256  ) ) 13  + (c156  + a (c56  + ) )u2>  u 

+ (C124  + “ (C24  + C224}  )U2,  22  + (C11  + C155  + C*C255)U3, 11  + (C144  + °rC244)U3,  221 
+ W1,2[C56U1,11+C24U1,22]  + W2, 1 (C14  + °56  )U2, 12  + E33  [ (C14  + C56  + C134 

+ c356  + P (C124  + C256 ) )U1, 12  + (C56  + °356  + P (C56  + °256 } )u2, 11 + (c24  + C234 
+ 3 (c24  + c224)  )u^  22  + (c13  + Pc12  + 2c55  + c355  + 3c255)u,^  1]l  + (2c44  + c344 

* Sc244>u3,22'  +*3,2I<'l2  + '«6,’,l,12+c66u2,U  + <’2j’,2,22  + 2cS6"3,U 
+ 2C24U3,22>  +"2,3tc44u2,22+c55“2,llI  * pS3  - 


where 


®11U1,11+  (®12 +e26)U2,12  + (®14 +e25)u3,12 +e26Ul,22 

<v 

Sll^,  11  “ ®22^22  = ° » 


or  = - ci2^C22  * ^ C-\VC')-)  • 


32  22 


(4.2) 


(4.3) 


From  (2.13),  with  (2.10),  for  straight-crested  waves  with  propagation  wave- 
number §,  we  have  the  boundary  conditions 

K2V  = 0 , ^2*"  ®o5*»  atx2  = 0,  (4.4) 

12  13 

and  from  (2.2)  and  (2.11)  we  obtain  y the  pertinent  constitutive  relations 

*21 = °56U3, 1 + C66  (ul,  2 + U2, 1>  + e26$,  2 + E11 1 (2c66  + C166  + *266  )ul,  2 
+ (c66  +C166  +Qf(c66  +c266))u2,l+  (c56  + C156  +“c256)u3,l1  +Wl,2t(c12 
+ C66)U1, 1 + C22U2,  2 + C24U3,  21  + W2, 1C66U2,  2 + E33  [ (C366  + *=266 )ul,  2 
+ (C366  + P(C66  +C266))U2,1  + (C56  + C356  +0C256)u3,l1  + W3,  2C66U3, 1 + W2,  3°56U2, 1 > 


*22  = C12U1, 1 + C22U2,  2 + C24U3,  2 + ei2?, 1 + E11 C (C12  + C112  + * (C12  + C122> )U1, 1 

+ (C122 +“(2C22 +C222))U2,  2 + (C124  + ® (°24  + C224)  * U3,  21  + "l,  2C22U1,  2 
+ W2, 1 f (C12  + C66 )U2, 1 + C66U1, 2 + C56U3, l1  + E33 [ (C123  + P (°12  + C122) )ul, 1 
+ (C223  + P(2C22+C222))U2,2+  (C24  + C234  + P(C24+C224))U3,2]  +W3,2C22U3,2 

+ W2,3[C14U1,1  + 2C24U2,2+C44U3,2J' 

*23  = C14U1, 1 + C24U2,  2 + C44U3,  2 + e14*  1 + Elll  (C14  + C114  + aC124)ul,  1 + (C124 

+ a(c24+C224))U2,2+  (C144+Q,C244)U3,2]  + ”l,  2C24U1, 2 + *2, 1C14U2, 1 

+ E33  [ (C14  + C134  + PC124)U1, 1 + (°24  ''  °234  + P (°24  + C224}  )U2,  2 + (2C44 

+ C344  + PC244)U3,23  + *3, 2 [C22U2, 2 + C12U1, 1 + 2C24U3, 21  + W2,3C44U2,2  ’ 

<\< 

\ = e25U3,l  + e26(U!,2+U2,l)  “ e22^,2  * (4’5) 

For  the  direct  calculations  performed  in  this  section  we  take  E*  and  E^3  to  be 
constant  with  respect  to  the  decay  length  of  the  surface  wave.  This  is  a 
reasonably  accurate  assumption  for  thick  substrates  for  which  the  surface 
wave  decays  rapidly  compared  to  the  variation  of  the  flexural  stress,  but 
not  for  relatively  thin  substrates.  In  addition,  in  the  direct  calculation 
performed  in  this  section  we  ignore  the  variable  displacement  gradient  terms 
W1  2’  W2  1 an<3  W3  2>  W2  3 in  and  (3*  14),  which  do  not  give  rise  to 

strains,  because  of  their  dependence  on  the  position  coordinate  X^.  The  in- 
fluence of  this  latter  assumption  is  not  too  great  because  the  neglected  terms 
do  not  contain  the  third  order  elastic  constants,  which  are  considerably  larger 
than  the  second  order  elastic  constants.  Calculations  are  performed  for 
cylindrical  flexure  either  in  the  X^-direction  or  in  the  X3~direction  only, 


not  in  both  directions  simultaneously.  Consequently,  in  the  calculations 


13. 


performed  the  biasing  terms  containing  either  the  index  3 or  1 in  (4.2) 
and  (4.5)  vanish. 

A solution  satisfying  (4,2)  and  (4.4),  with  (4.5),  may  be  taken  in  the 
usual  form 


- I 


„«,<»)  ie«?X2 
C A.,  e e 


m=l 


4 

5- 1 

m=l 


, . a»i  i§OC.-Vt) 

c(m)B(m)e  m 2fi  1 


(4.6) 


where  V,  3 , A^m\  B ^ and  C ^ are  determined  numerically  so  that  (4.2) 
nr  V 7 

14-16 

and  (4.4),  with  (4.5),  are  satisfied  in  the  usual  manner  . Calculations 

17  18 

have  been  performed  for  ST-cut  and  Y-cut  quartz  7 in  X^-flexure  and  the 

results  are  plotted  as  the  line  labeled  a in  Figs. 2 and  3,  respectively.  In 

addition,  we  have  performed  calculations  for  ST-cut  quartz  in  X^-f lexure  when 

the  biasing  terms,  i.e.,  those  proportional  to  E^,  have  been  ignored  in  the 

boundary  conditions  (4.4),  with  (4.5),  and  the  resulting  line  is  labeled  a* 

in  Fig. 2.  it  is  clear  from  the  figure  that  the  difference  is  quite  great  and 

that,  ‘consequently,  in  general  the  nonlinearities  in  the  boundary  conditions 

are  quite  important.  Calculations  have  been  performed  for  X^-flexure  also, 

and  the  results  appear  as  the  line  labeled  b in  Figs. 2 and  3 for  ST-cut  and 

Y-cut  quartz,  respectively.  No  line  labeled  a*  appears  in  Fig. 3. 

For  the  case  of  static  flexural  biasing  stresses  both  in  and  normal  to 

the  propagation  (X3) -direction  for  Y-cut  lithium  niobate,  Eqs.  (2.1)  and  (2.2), 

12  13 

with  (2.11),  the  appropriate  relations  in  Sec. 3 and  (4.1)  take  the  form  ’ 


14 


C22U2,  22  + (C23  + C44)U3,  23  + 2C24U2,  23  + C24U3,  22  + C44U2,  33  + *22$,  22^  (C32  + e24)CP,  53 
+ E33 1 (C223  + P (2C22  + °222) )U2, 22  + (2°234  + P (4°24  + 2C2245 )U2, 23  + <C33  + °344 
+ P (2c44  + C244} >U2, 33  + (C24  + C234  + P (°24  + C224) ,U3, 22  + (C23  + °44  + °233  + °344 
+ P (C23  + C44  + °223  + C244) )U3, 23  + (°34  + C334  + P (C34  + C234) >U3, 331 
+ W3, 2 [C22U3, 22  + 2C24U3, 23  + °44U3, 331  + W2, 3 t2C24U2, 22  + 2 (C23  + C44)U2,  23 

+ 2c34U2,33  +C44U3,  22  + 2°34U3,  23  + C33U3,33]  + E11 [ (C122  + C*  (2C22  + C222}  )U2,  22 

+ {2ci24  + a(4C24  + 2C224))U2,23+  (C13  + C144  + “ (C23  + 2°44  + W )U2,  33  + (C124 
+ ® (C24  + C224) )U3, 22  + (C123  + C144  + “ (C23  + C44  + C223  + C244 > )U3, 23  + (C134 

+a(C34  + C234))U3,331  = pU2 » 

/v> 

C24U2, 22  + (C23  + °44)U2, 23  + °44U3,  22  + °33U3, 33  + e24^, 22  + e33^, 33  + E33 [ (°24 


+ C234  + P (C24  + C224) )u2, 22  + (°23  + C44  + °233  + C344  + P (°23  + °44  + C223 

+ C244} )U2, 23  + (C34  + °334  + P (°34  + C234) >U2, 33  + (2C44  + °344  + P°244)U3, 22 

+ (4C34  + 2C334  + 2PC234)U3, 23  + (2°33  + °333  + °33  + Pc233)U3, 33  + *3, 2 tc22U2, 22 

+ 2C24U2, 23  + °44U2, 33  + 2°24U3, 22  + 2 (°23  + °44)U3, 23  + 2C34U3, 33J  + W2, 3 [C44U2, 22 

+ 2C34U2,23+C33U2,33]  + En  [ (C124  + ® (C24  + C224  > )U2,  22  + (C123  + C144  + " (C23 
+ C44  + C223  + W ' )U2,  23  + (C134  +“  <C34  + C234))U2,  23  + (C144  + “C244)  U3,  22 

+ 2(C134  + QC234)U3,23+  (c13  + aC23  + C133  + QC233)u3,33]  = P*3  ’ 

e22U2, 22  + (e24  + 632)u2, 23  + 624U3,  22  + 633U3, 33  " *22*  22  " 833* 33  = ° > (4‘ 7) 


|l 


where  O'  and  |3  are  given  in  (4.3).  Again  the  boundary  conditions  are  given  by 

12  13 

(4.4)  and  from  (2.2)  and  (2.11)  we  obtain  ’ the  pertinent  constitutive 
equations 

*22  = C22U2, 2 + C23U3, 3 + °24  (U2, 3 + U3,  2 * + ®22^,  2 + ®32^, 3 + E33  [ (°223  + ^ (2°22 
+ °222) )U2, 2 + (C234  + ^ (2°24  + °224) )U2, 3 + (°24  + C234  + ^ (c24  + °224) )u3, 2 
+ (C23  + C233  + P (C23  + C223) )u3, 31  + W3, 2 tc22U3, 2 + °24U3, 31  + W2, 3 [2c24U2, 2 
+ (C23+C44)U2,3+C44U3,2+C34U3,3]  + E11 [ (C122 + “ (2°22 + C222> >U2,  2 


+ (c124  + Of  (2c24  t c224)  )u2j  3 + (c124  +a(c24  + c224)  )u3>  2 + (c123 

+ a(C23+C223))U3,3]  > 

*23  = C24U2, 2 + °34U3, 3 + °44  (U2,  3 + U3,  2 ’ + ®24$, 2 + E33 1 (C24  + C234  + P (C24 
+ C224 > >U2, 2 + (°44  + C344  + P (C44  + C244) )U2, 3 + (2c44  + C344  + eC244)U3, 2 
+ <2C34  + C334  * PC234)U3, 31  + ” 3 , 2 fC22U2,  2 + C24U2, 3 + 2C24U3, 2 + (C23 
+ C44)U3, 31  + W2, 3 [C34U2, 3 + C44U2, 21  + EJl 1 1 (C124  + “ (C24  + C224) )U2, 2 


+ (C144+a(C44+C244))U2.3+  (C144+“C244)U3.2+  (C134+<>C234)U3.31  » 


°2  ®22U2,2 +e24(U2,3  + U3,2)  " S22Cp,  2~  e23Cp,  3 ‘ 


A solution  satisfying  (4.7)  and  (4.4),  with  (4.8),  may  be  taken  in  the 


usual  form 


UV  = L C \ e e 

m=l 


, V * 1 


i§(x,-vt) 

c(m)B(m)e  m 2 e 3 


r 


where  V,  3^,  and  C ^ are  to  be  determined  numerically  so  that  (4.7) 

15  18  19 

and  (4.4),  with  (4.8)  are  satisfied  in  the  usual  manner  ’ ’ . Calculations 

20,  7 

have  been  performed  for  Z-propagating  surface  waves  on  Y-cut  lithium  mobate  , 
and  the  resulting  lines  are  plotted  in  Fig. 4.  The  notation  on  the  lines  in 
Fig. 4 has  the  same  meaning  as  in  the  case  of  quartz  in  Figs. 2 and  3,  and  no 
a*  line  is  plotted. 


5.  Perturbation  Procedure 

If  we  define  the  change  in  phase  velocity  due  to  the  applied  flexural 


biasing  state  by 


AV  = V - v 


then  it  has  been  shown  that  for  traction  free  and  zero  normal  component  of 
electric  displacement  boundary  conditions  the  first  perturbation  in  phase 


velocity  is  given  by' 


av  " V2v«5m  « 


where 


® 2TT/S, 

H„~J  J "C"",!.  * <1X2dXl 


~ A AM  , A ~M 

KLY_ClW\>,K  eKL-^,K  > 

A A 

'vn  a aM  a a/M 

~ eLK*Ua,  K ” eLKCp,  K ’ 


«M  M . ~M  ~M  . 

UY  = V V =tf>  /NM  ' 


2 of  f MM 

NM  = P J J UVUY  *2*1  ’ 


where  the  coefficients  in  (5.4)  are  defined  in  (2.3).  In  the  case  of  a purely 
elastic  bias,  which  exists  when  electric  and  electroelastic  coefficients  are 


17. 


ignored,  (2.11)  and  (2.12)  hold  in  place  of  (2.3),  and  in  place  of  (5.3)  and 
(5.4)  we  have 

00  2rr/£ 

v- J 


(5.6) 


1 * * M 

=*  c u 

jY  LYKV  v,  k ’ 


(5.7) 


where  is  given  in  (2.11). 

When  (5.6)  is  applied  to  the  case  of  rotated  Y-cut  quartz  subject  to 

cylindrical  flexure  in  the  propagation  (X] ) -direction,  we  obtain  the  form 
* 2tt/5m 

«„"K  J r*Liai,  1 + ^252, 1 * %353, 1 + *2l”l,  2 


+ ^2G2,2  +^3a3,21 


(5.8) 


The  substitution  of  (5.7),  with  k2.11),  along  with  the  appropriate  relations 
in  Sec. 3 into  (5.8)  yields 

2tt/A 


IV-JdX2  J dXl 


[ (cIl  + 2C11  + Clll  + ac112)Eu\  l\  1 + (CIl  + C166 


+ a(2c66+C266))Ell{i2,l(i2,l+  (cU  +°155  + ac255,El1G3>  l\l  + (2°66 
+ °166 +0(C266)E11U1,  2U1,2  + (°122  + “ (2°22  + C2225  )E11U2,  2U2,  2 + (°144 
+ OC244)EllU3, 2U3, 2 + 2 (c12  + C112  + “ (cl2 + °122) )E11Q1, 1^2, 2 + 2 (°14 


+ C114+“C124)E11U1,1U3,2  + 2(C66  + °166  + “ (c66  + c266 } )Ellul,  2U2, 1 
+ 2,C56  + c156  t<'c256,Eil5l,2il3,l  + 2(<:lS6  + “ lcS6  + C256>  >E}l53,  A,  1 
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where  from  (3.10)  and  (3.11),  with  Fig.l 


_(i>  a (i  _ h.)  es  „ . !!i  Es 

Eii  V-  TV  E ' w2,i  ir E » \ 


x, 

-i* 


(5.10) 


ES«-  3M/2h2c*1 


(5.11) 


When  the  integrals  in  (5.9)  are  performed,  the  resulting  expression  is 

extremely  cumbersome  and  not  terribly  revealing.  Consequently,  we  do  not 

present  the  complete  equation  from  which  the  calculations  are  made.  However, 

in  order  to  give  some  idea  of  the  type  of  terms  that  occur,  we  present  two 

typical  distinct  types  of  integrals.  These  are 
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the  * denotes  complex  conjugate  and  2h  is  the  thickness  of  the  plate,  which  is 


many  times  larger  than  the  wavelength  of  the  surface  wave  so  that 


(u^,  cp)  -*  0 at  Xj  = 2h  . (5.15) 

Note  that  the  biasing  term  in  the  typical  integral  in  (5.13)  is  the  type  of 
term  that  was  ignored  in  the  direct  calculations  performed  in  Sec. 4.  When  all 
terms  of  this  type,  i.e.,  those  not  giving  rise  to  a strain,  are  ignored  along 
with  the  X2~dependence  of  the  strain  in  the  perturbation  calculations  for 
both  ST-cut  and  Y-cut  quartz,  the  results  are  indistinguishable  from  the  lines 
labeled  a in  Figs. 2 and  3,  respectively.  However,  when  the  spatial  variation  of 
all  biasing  terms  in  (5.9)  is  included  in  the  perturbation  calculations  for  ST- 
cut  and  Y-cut  quartz,  the  resulting  curves  plot  as  the  lines  labeled  a7  in  Figs. 2 
and  3,  respectively.  It  can  be  seen  from  the  figures  that  although  the  difference 
between  curves  a and  a7  is  small,  for  relatively  thick  substrates,  it  increases 
appreciably  with  decreasing  thickness-to-wavelength  ratio.  This  means  that  since 
the  change  in  such  physical  quantities  as  phase  velocity  with  biasing  deformation 
is  very  small,  the  perturbation  procedure  is  significantly  more  accurate  than  the 
direct  calculation  for  the  determination  of  such  quantities  because  it  enables 
the  influence  of  variable  biasing  terms  to  be  included  in  the  calculation. 

When  perturbation  calculations  are  performed  for  ST-cut  and  Y-cut  quartz  subject 
to  cylindrical  flexure  in  the  X^-direction,  i.e.,  normal  to  the  propagation  direc- 
tion, the  resulting  curvet'  plot  as  the  curves  labeled  b7  in  Figs.  2 and  3,  respect- 
ively. In  this  case  no  terms  due  to  the  variable  biasing  displacement  gradients 
w2  3 and  w3  2 survive  in  the  perturbation  calculation  because  they  are  odd  in  X3 
and  the  surface  wave  is  independent  of  X3> 


When  (5.6)  is  applied  to  the  case  of  Y-cut  lithium  niobate  subject  to 
cylindrical  flexure  in  the  propagation  (X3 ) -direction,  we  obtain  the  form 
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The  substitution  of  (5.7),  with  (2.11),  along  with  the  appropriate  relations 
in  Sec. 3 into  (5.16)  yields 

7 "tfM  . 

HM-“  J **2  J <«3[(C223  + 0(2C22+C222))E33ii2,2"2f2  + 2(C234  + P(2C24 

O 

+ c224)  >=33*2,  2*2,  3 + 2 (C.24  + C234  + P (c24  + C224)  )E33G2,  A,  2 + 2 (C23  + C233 
+ P (c23  + c223) )E33u^  2u3^  3 + (c33  + c344  + P (2c44  + c244) )E33u2^  3u2^  3 + 2 (c44 

+ C344  + P (C44  + C244) )E33U2, 3U3, 2 + 2 (C334  + Pc234)E33U2, 3U3, 3 + (2°44  + °344 
+ eC244)E33U3, 2U3, 2 + 2 (C334  + 0C234)E33U3, 2U3, 3 + (C33  + 2°33  + C333 
+ ^C233)E33U3, 3U3, 3 " 2C24W3, 2U2,  2U2, 2 “ 2 (C23  + °44)W3, 2U2, 2U2, 3 + 2 <C22 
- C44)W3, 2S2, 2*3, 2 + 2c24W3, 2*2,  2*3,  3 + 2°24W3,  2*2,  3*3,  2 + 2 (c44 


" C33)W3, 2°2, 3U3,  3 + 2°24W3, 2U3,  2U3,  2 + 2 (C23  + C44)W3,  2U3,  2U3,  31  ’ t5- 17) 


where  from  (3.14)  and  (3.15),  with  Fig.l,  we  have 
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As  in  the  case  of  quartz  and  for  the  same  reasons  we  do  not  present  the  full 
equations  from  which  the  calculations  were  made,  but  instead  present  the  two 


typical  distinct  types  of  integrals 
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where 
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The  results  of  perturbation  calculations  for  Y-cut,  Z-propagating  surface  waves 
yield  curves  a;  and  b'  in  Fig. 4. 
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FIGURE  CAPTIONS 


Figure  1 Schematic  Diagram  of  the  Surface  Wave  Structure  Showing 

the  Applied  Static  Loading 

Figure  2 Relative  Change  in  Surface  Wave  Velocity  Per  Unit 

Applied  Maximum Extensional  Biasing  Strain  vs.  Thickness- 
Wavelength  Ratio  for  ST-X  Quartz.  The  curves  a'  and  b' 
are  for  cylindrical  flexure  in  and  normal  to  the 
direction  of  propagation,  respectively. 

Figure  3 Relative  Change  in  Surface  Wave  Velocity  Per  Unit 

Applied  Maximum  Fxtensional  Biasing  Strain  vs.  Thickness 
Wavelength  Ratio  for  Y-X  Quartz.  The  notation  is  the 
Same  as  in  Fig. 2. 

Relative  Change  in  Surface  Wave  Velocity  Per  Unit 
Applied  Maximum  Extensional  Biasing  Strain  vs.  Thickness 
Wavelength  Ratio  for  Y-Z  Lithium  Niobate.  The  notation 
is  the  same  as  in  Fig. 2. 


Figure  4 


